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Singularities in the dark energy late universe are discussed, under 
the assumption that the Lagrangian contains the Einstein term R 
plus a modified gravity term R a , where a is a constant. The 4D 
fluid is taken to be viscous and composed of two components, one 
Einstein component where the bulk viscosity is proportional to the 
scalar expansion 9, and another modified component where the bulk 
viscosity is proportional to the power 8 2a ~ 1 . Under these conditions it 
is known from earlier that the bulk viscosity can drive the fluid from 
the quintessence region (w > —1) into the phantom region (w < —1), 
where w is the thermodynamical parameter [I. Brevik, Gen. Rel. 
Grav. 38, 1317 (2006)]. We combine this 4D theory with the 5D 
Randall-Sundrum II theory in which there is a single spatially flat 
brane situated at y = 0. We find that the Big Rip singularity, which 
occurs in 4D theory if a > 1/2, carries over to the 5D metric in the 
bulk, \y\ > 0. The present investigation generalizes that of an earlier 
paper [I. Brevik, larXiv:0 807.1797: to appear in Eur. Phys. J. C] in 
which only a one-component modified fluid was present. 

Keywords: modified gravity; viscous cosmology; Randall-Sundrum model 

1 Introduction 

There has recently been a great deal of interest on modified gravity theories 
in four dimensions, and what consequences they may have for the Big Rip 
singularity in the late universe (cf., for instance, the review of Copeland 
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et al. [1J). The equation of state for the cosmic fluid is conventionally 
written as p = wp = (7 — l)p, where w = — 1 corresponds to a vacuum 
fluid (or cosmological constant), -1 < w < -1/3 to a quintessence fluid, 
and w < — 1 to a phantom fluid. The last-mentioned case predicts a Big 
Rip singularity in the late universe. Both quintessence and phantom fluids 
lead to the thermodynamical inequality p + 3p < 0, thus breaking the strong 
energy condition. 

It is of interest to combine the modified gravity in 4D with the Randall- 
Sundrum model in 5D. We shall consider the RS II model, in which there 
is a single spatially flat brane situated at y — 0, surrounded by an AdS 
space [2J. Therewith one can analyze the relationship between the 4D Big 
Rip singularity and the corresponding behavior of the metric in the 5D bulk, 
\y\ > 0. This task was recently undertaken in Ref. [3], under the assumption 
that there was a one-component uniform modified viscous fluid on the brane. 
The bulk viscosity ( was taken to vary with the scalar expansion 6 as Q oc 
6* 2a_1 , with a is a constant present in the 4D action 



Here k\ = SnG^, f is a constant, and L m is the Lagrangian of the matter. As 
shown in Ref. [I] , when ( varies in this way the formalism leads naturally to 
a Big Rip singularity in a finite future time. (This kind of theory is a natural 
generalization of the theory for a viscous Einstein fluid, a = 1, as delineated 
in Ref. [5].) The main result of the investigation in Ref. [3] was that the 4D 
singularity on the brane becomes carried over to the 5D bulk, \y\ > 0. The 
scale factors on the brane and in the bulk are actually proportional to each 
other. 

The new element in the present investigation is that we consider a two- 
fluid model, with action 



It means that we include an Einstein fluid in addition to the modified fluid. 
There are physical motivations for this generalization. As emphasized by 
Vikman, for instance, [6J, it is necessary to introduce a two-component model 
within the framework of the scalar field picture. There are moreover several 
recent works on cosmological models with two scalar fields; one may consult 




(1) 
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Refs. |3 El El EH m US [131, for instance. Equation © is better motivated 
physically than Eq. ([!]). 

One may ask: is there the same close relationship in the two-fluid case 
between the singularities on the brane and in the bulk as in the single-fluid 
case? The answer turns out to be affirmative, as we shall show below. 

Finally, we mention that there exist more general f(R) theories. A general 
review of modified gravity can be found in Ref. [14], and recent reviews of 
f{R) gravity theories unifying dark energy, inflation, and dark matter, can 
be found in Refs. [T5l [16] .) 

2 A brief resume of the 4D theory 

We begin by reproducing some salient features of the 4D theory for a two- 
fluid system, following the treatment in Ref. [17] . The spatially flat FRW 
metric is 

ds 2 = -dt 2 + a 2 (t)dx 2 , (3) 
and the energy-momentum tensor of the viscous fluid is 

T,„, = pUyUv + ph^ u , (4) 

where h^ v = g^ v + U^U U is the projection tensor and p = p — (9 the effective 
pressure. In comoving coordinates, U° = 1, U l = 0. From variation of the 
two-fluid action (T5]) we can derive the equations of motion (we put A4 = 0). 
Of main interest is the (OO)-component of the equations. We combine that 
particular component with the energy conservation law 

p + (p + p)3H = %H 2 , (5) 

which in turn is a consequence of the covariant conservation equation V U T^ U = 
0. Here H = a/a = 9/3 is the Hubble parameter. Some calculation leads to 
the equation 

6H + 9 7 # 2 + ^f R a ~ 3a/ [(37 - 2)H + 3 1 H 2 ]R a ~ 1 

+ 3a(a - l)/ [(37 - l)HR + R]R a ~ 2 + 3a(a - l)(a - 2)f R 2 R a ~ 3 = 9kI(H, 

(6) 

which is a nonlinear differential equation for H(t) in view of the relationship 
R = Q(H + 2H 2 ). 
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Equation is complicated. We will look for solutions of the form 

H = —t, where X = 1 - BH* t. (7) 
A 

Here B is a nondimensional constant, whose value has to be calculated on 
the basis of initial assumptions for w, a, and ( for the two fluid components. 
The star subscript refers to the initial (present) time £* = 0. If B > 0, H 
will diverge in a finite singularity time t = t s , and Big Rip will occur. 

Assume now that the total bulk viscosity ( is a sum of two parts, one 
part (e referring to the Einstein fluid and a second part ( a referring to the 
modified fluid. As mentioned above, they will be taken to be proportional 
to 9 and 9 2a ~ 1 , respectively. The corresponding proportionality factors are 
denoted te and r a . Thus 

Ce = St e H j ( a = T a (3H) 2a -\ (8) 

A consequence of these assumptions is that Eq. ([6]) becomes satisfied for 
the Einstein component and the modified component separately; the time- 
dependent terms automatically drop out. It turns out that there is a re- 
lationship between the factors t e and r a , the form of r a = T a {r E ) being 
determined from a and /q as well as from w. The energy conservation law 
(j3J) holds for each component separately. From these equations we get two 
different expressions for the constant B: 

B = -?i + ^ (9) 

* * P*E 

2a 2a p* a 

We have here for t = the relations p* = p*E + P*a', moreover 

3H*\ 2a 1 p^E P*a /,,n 

> PE=x2, P*=X^- ( U ) 

3 Implications for the 5D theory 

Now move on to consider the RS flat brane situated at y — 0, surrounded by 
an AdS space with metric 

ds 2 = -n 2 (t, y)dt 2 + a 2 (t, y)5 ij dx i dx j + dy 2 . (12) 
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Here n(t, y) and a(t, y) are determined from Einstein's equations 



1 2 
Rab — ^9abR + gAB^ = k Tab, 



(13) 



where k 2 = 87TG5 . Note that whereas we put the 4D cosmological constant 
equal to zero above, we keep the 5D cosmological constant A(< 0) different 
from zero in order to comply with the main idea of the RS model. The 
coordinates are x A = (t, x 1 , x 2 , x 3 , y). From the junction conditions across 
the brane we get for arbitrary y, after integration with respect to y [3], 



a 

na 



c_ 

7a' 



(14) 



On the brane we may take n Q (t) = 1. Omitting the unimportant C term (the 
"radiation term"), we obtain then on the brane 



I K 4 

-AH (a + p) 2 , 

6 36 V H! 



(15) 



where a is the brane tension, assumed constant. We let generally the sub- 
script zero refer to the brane. 

The 5D equation ( JT5l) is important in our context, as it may be regarded to 
represent a bridge between the 4D and 5D cosmologies. Namely, by inserting 
for p = Pe + Pa from Eq. (TTTj) . we get 



Hi 



1 , K 

-Ah — 

6 36 



P*E 



be, ty 



+ 



BH,t) 



2a 



(16) 



We shall investigate the behavior of the scale factor near the Big Rip, occur- 
ring at t s = 1/ (BH*). We then assume that B is a positive quantity. We get 
approximatively 



Oo 

a 



K 



P*E 



+ 



P* 



l-BH*t) 2 (l-BH*t) 2a 



(17) 



showing that the constants A and a are no longer important in this limit. 
The solution is of the form 



a (t) ~ exp 



{K 2 /6)p* 1 



(BH*) 2 (t s -t) (2a-l)(BH^(t s -t) 2 ^ 



(18) 



5 



Two different sub-classes need here to be distinguished: 



(%) The case a < 1. The influence from the modified fluid component then 
becomes subdominant near the Big Rip. The behavior of a(t) is governed by 
the Einstein fluid component. It might seem natural to suggest that this is 
after all the most likely scenario in the late universe. 

(ii) The case a > 1. Then, the modified fluid component will dominate 
at the end, regardless of what was the initial ratio between p* E and p* a at 
the initial instant t = 0. Moreover, the strength of the singularity is seen to 
be larger than in the Einstein case, the strength increasing for larger input 
values of a. 

Let us now go back to the 5D equation (fT4l) in the bulk. Taking into 
account the relation 

n(t,y) = ^, (19) 
a (t) 

which in turn is a consequence of there being no energy flux in the y direction 
from the brane T ty = (cf., for instance, Ref. [18]), we obtain the solution 



(20) 

with p = ^J—A/6. (Recall that we have assumed C = 0, k = 0.) This is 
formally the same solution as in the single- fluid case [3]. The characteristic 
properties of the two- fluid system turns up in the prefactor ao(t), but they 
do not influence the variation of a(t, y) upon y in the bulk. 

We may thus summarize our work as follows: 

• The Big Rip singularity on the brane as following from the 4D theory 
- regardless of whether it is the Einstein fluid or the modified fluid that 
is the dominant component - becomes immediately transferred to the bulk, 
\y\ > 0. The brane tension a and the 5D cosmological constant A play no 
role for the establishment of the brane singularity, but they turn up again 
in the dependence of a(t,y) upon y in the bulk region; cf. Eq. (120]) . The 
reason for this immediate effect is not known, but may be related to quantum 
mechanics. 
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• Our basic action (j2J) with a constant means a simple example of modified 
gravity. We assumed here a straightforward superposition of an Einstein fluid 
and a modified fluid. (More complicated modified gravity theories can be 
found, for example, in Refs. [HI IS, US].) An important point in our analysis 
was to assume both fluid components viscous, and to assume the specific 
forms given in Eq. (jSJ) for the bulk viscosities. These specific forms admit 
both fluid components to pass through the w = — 1 barrier into the phantom 
region, and thereafter into the Big Rip singularity. 

• The nature of the singularity on the brane was shown to depend on 
the magnitude of a. If a < 1, the singularity near Big Rip was determined 
by the Einstein component, whereas for a > 1 it was determined by the 
modified component. This kind of behavior was independent of the relative 
magnitude of the energy densities p*E and p* a at the initial instant t — 0. 
When the modified component dominates, the singularity is always stronger 
than in the Einstein case, and the singularity becomes stronger the higher is 
the value of a. As simple examples one may note that the value a = 1/2, 
often considered, corresponding to a y/~R term in the Lagrangian, belongs to 
the Einstein case whereas the value a = 2 (a R 2 term in the Lagrangian) 
belongs to the modified case. 

• It should be emphasized that we have assumed the bulk viscosity only, 
simply omitting the shear viscosity which is known to be the most important 
component in ordinary hydrodynamics. This seems to be a natural way to 
proceed, all the time that we are assuming spatial isotropy in the cosmic 
fluid. 
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